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Section  1 
Introduction 

Th  is  paper  describes  the  theoretical  foundations,  development,  and  verification  of  two 
discrete-element  processors,  SWEEPS  and  TARGET,  that  can  treat  the  steady  state 
acoustic  radiation  and  scattering  of  a  resilient  submerged  structure.  Both  are  built 
around  a  database  associated  with  the  Underwater  Shock  Analysis  (USA)  code  [l]  as 
they  share  a  common  approach  to  the  fluid-structure  interaction:  the  use  of  a  boundary- 
element  method  based  on  the  Doubly  Asymptotic  Approximation  (DA A)  [2]. 

DA  A  methods  have  long  proved  their  usefulness  in  underwater  shock  problems  and  it  is 
only  recently  that  they  have  been  applied  to  acoustics.  The  use  of  several  forms  of  the 
D.^A  to  study  the  forced  vibration  of  a  submerged  spherical  shell  [3]  has  shown  that 
the  first  order  DAAi  is  generally  unsatisfactory  for  acoustics  applications,  however  tw'o 
second  order  forms  appear  to  perform  quite  well  in  most  cases.  These  are  the  so-called 
curvature  corrected  form  DAAj  and  the  modal  form  DAA^  . 

Although  theoretically  exact  solutions  to  the  underwater  acoustics  problem  do  exist 
\4-7],  the  DA.\  approach  uses  a  fluid  mass  matrix  and  a  diagonal  fluid  area  matrix  that 
do  not  depend  upon  frequency.  This  fact  leads  to  a  more  efficient  use  of  computational 
resources  when  performing  variable  frequency  calculations.  In  contrast,  the  governing 
fluid  matrices  must  be  reformed  in  full  for  every  frequency  in  the  exact  formulations.  In 
addition,  the  exact  solution  can  be  prone  to  the  well  known  critical  frequency  problem 
;4,6]  that  the  DAA  approach  does  not  encounter,  although  the  exact  method  discussed 
in  [7]  does  avoid  the  problem.  Hence,  it  appears  that  DAA  methods  can  be  used  to 
advantage  in  underwater  acoustics  in  that  an  increase  in  efficiency  can  offset  some  loss 
of  accuracy.  Indeed,  if  preliminary  results  can  be  obtained  with  the  use  of  D.\.42  or 
DAAJ*  and  some  frequency  ranges  prove  to  be  of  particular  interest,  then  an  exact 
approach  could  always  be  used  in  this  limited  range. 

The  primary  aim  of  this  paper  is  to  document  the  problem  formulation  used  in  the  new 
processors  and,  to  verify  the  integrity  of  the  software  by  solving  both  a  radiation  and 
a  scattering  problem  for  the  spherical  shell.  For  information  on  the  operation  of  the 
code  itself,  and  sample  input  and  output,  a  usage  primer  has  been  written  [8). 

In  the  next  section  the  governing  finite-element,  boundary-element  equations  are  pre¬ 
sented  for  the  steady  state  vibration  of  a  submerged  structure  excited  either  by  a  set  of 
internal  forces  with  the  same  specified  frequency  but  otherwise  arbitrary  magnitudes 
and  phases,  or,  by  an  infinite  train  of  sinusoidal  incident  waves  emanating  from  a  spher¬ 
ical  source  with  a  specified  frequency  and  magnitude.  A  simple  and  direct  elimination 
solution  is  then  given  for  the  structural  displacement  field  and  the  wet  surface  scattered 
pressures.  This  solution  process  forms  the  basis  of  the  SWEEPS  processor. 

The  following  section  briefly  discusses  the  use  of  the  Helmholtz  integral  equation  in 
order  to  find  fluid  pressures  away  from  the  wet  surface  of  the  structure.  A  discretized 
form  of  this  equation  is  implemented  in  the  TARGET  processor. 
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The  next  section  summarizes  the  governing  equation  system  and  selected  analytical 
DAA2  solutions  for  the  spherical  shell.  These  are  then  called  upon  in  the  code  verifi¬ 
cation  section  which  follows.  The  discrete-element  computational  results  are  found  to 
exhibit  excellent  agreement  with  the  closed  form  solutions.  Finally,  some  recommen¬ 
dations  for  future  work  are  presented  in  the  last  section. 
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Section  2 

Governing  Equations  For  Wet  Surface  Solution 


The  interaction  equations  for  a  DAA2  time-harmonic  vibration  analysis  of  a  submerged, 
linear-elastic  structure  may  be  written  in  matrix  form  as  |2] 


where 


E^ii  —  — +  iujCs  +  Kfi, 

Es/  =  G  A/, 

E/a  =  pciuj^  (wM /  -  ifl /M /)G^, 
E//  =  — f  +  pc  (icjA /  -|-  n /A /), 

6a  ~  ~ 

gf  =  pcu!  (wM /  -  7X1  /M /)u;. 


(1) 


(2) 


Here  Ma,  Ca,  and  Ka  are  the  structural  mass,  damping,  and  stiffness  matrices,  respec¬ 
tively,  G  is  the  fluid-structure  transformation  matrix  relating  fluid  node  point  forces 
normal  to  the  wet  surface  of  the  structure  to  node  point  forces  in  the  structural  compu¬ 
tational  system,  A /  is  the  wet  surface  fluid  element  area  matrix,  M /  is  the  wet  surface 
fluid  mass  matrix,  and  fl/  is  the  wet  surface  fluid  frequency  matrix;  x  is  the  structural 
displacement  vector.  P5  is  the  wet  surface  scattered  pressure  vector,  is  an  applied 
structural  force  vector,  py  is  an  incident  wave  w’et  surface  pressure  vector,  and  U/  is 
an  incident  wave  wet  surface  normal  velocity  vector.  (If  py  -  uy  -  0,  the  scattered 
pressure  P5  reduces  to  the  radiated  pressure  Pyy.)  In  addition,  p  and  c  are  the  fluid 
density  and  sound  speed,  respectively,  i  =  and  u>  is  the  frequency  of  steady  state 

vibration.  A  superscript  T  denotes  matrix  transposition. 

The  real,  symmetric  matrices  M^,.  C^,  and  Kj,  can  easily  be  generated  by  any  finite- 
element  structural  analysis  code,  and,  in  the  work  reported  upon  here,  STAGS  (STress 
Analysis  of  General  Shells)  [9^  has  been  used.  The  real,  diagonal  matrix  A/  is  trivially 
obtained,  while  the  real,  symmetric  matrix  M/  can  be  computed  by  the  boundary- 
element  method  of  [lOj.  These  two  fluid  arrays  ais  well  as  the  real,  rectangular  trans¬ 
formation  matrix  G  are  already  produced  by  the  FLUMAS  processor  of  the  USA  code. 
Finally,  the  real  matrix  fl/  may  be  obtained  from  either  of  two  formulations,  and  is 
such  that  the  matrix  product  fl/A/  is  symmetric.  The  development  of  [2],  which  is 
based  upon  the  method  of  fluid  boundary  modes,  gives 

nj  =  gpcAfMj\  (3) 
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where  </  is  a  scalar  parameter  that  can  vary  between  zero  and  unity,  g  =  0  reduces  (2) 
to  the  DAAi  equations,  g  =  1/2  appears  to  be  best  for  the  infinite  cylindrical  shell, 
and  g  =  1  is  best  for  the  spherical  shell.  On  the  other  hand,  the  formulation  of  Jll) 
does  not  contain  any  arbitrary  parameters  as  in  (3).  It  is  based  upon  the  method  of 
matched  asymptotic  expansions,  and,  for  the  fitting  procedure  described  in  [3],  yields 

n}  =  pcAjMJ^  ~cK,  (4) 

where  1C  is  a.  diagonal  matrix  of  wet  surface  mean  local  curvatures.  It  should  be  noted 
that  both  (3)  and  (4)  do  not  involve  any  additional  information  that  is  not  already 
provided  by  the  FLUMAS  processor,  in  particular,  the  mean  local  curvatures  are  used 
in  the  computation  of  M/  10]. 

For  convenience,  (3)  and  (4)  can  be  combined  into  the  one  expression 

flf  =  gpcAfMJ^-ScK,  (5) 

where,  if  ^  -  0  the  DAA^  form  is  obtained;  whereas,  if  6  =  g  =  1  the  DAAj  form  is 
obtained.  With  this  substitution,  (2)  then  become 

Eas  =  — U/'Ma  +  iujCf  -r-  Ka? 

Ea/  =  G  A/, 

Eifs  =  pciu)"^  \{u  +  iScK)^./ -  igpcAf]G^ , 

Ey/  =  -w^M j  -r  pc  -  ScK)A  f  -  gpcD /], 
ga  =  fs  -  G  A/Py, 

gy  =  pcijj  [(tj  +  i^cK')M y  -  igpcA y  Uy, 
w’here  Dy  is  the  symmetric  matrix  given  by 

Dy=AyMJ^A.  (7) 

Perhaps  the  most  important  characteristic  of  (6)  is  that  the  matrices  Ma,  Ca,  Ka,  G, 
Ay,  My.  an  d  Dy  are  frequency-independent,  so  that  they  need  only  be  computed  once 
for  a  complete  set  of  frequency-sweep  calculations.  This  characteristic  also  renders  (l) 
particularly  amenable  to  incremental  iterative  methods  of  solution  thus  avoiding  costly 
refactorization  of  the  coefficient  matrices  at  every  frequency  step.  More  will  be  said 
about  this  later. 

To  complete  the  governing  equation  system,  the  form  of  the  right  hand  side  forcing 
vectors  in  (2)  must  now  be  specified.  The  elements  of  the  internal  forcing  vector  fg  can 
be  written  ais 

=  (8) 
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where  F,  and  Oi  are  the  magnitude  and  phase  angle  respectively  of  the  degree  of 
freedom  of  fg.  Also,  the  elements  of  the  external  forcing  vectors  py  and  U/  can  be  given 
for  a  train  of  spherical  incident  waves  as 


Pii=Po^e  '  ‘  f. 

u/,  ^  —’(1  -  t/kRi)')i. 

pc 


Here  S  is  the  standoj  .  i.e.,  the  distance  between  the  origin  of  the  spherical  wave  and 
the  nearest  point  on  the  wet  surface  of  the  structure,  i?,  is  the  distance  from  the  origin 
of  the  spherical  wave  to  the  fluid  node  on  the  wet  surface,  and,  'y,  is  the  cosine  of 
the  angle  between  the  vector  corresponding  to  /?,  and  the  wet  surface  outward  normal 
vector  at  the  fluid  node,  po  is  the  amplitude  of  the  incident  pressure  at  the  standoff 
distance,  and  k  is  the  wave  number  w/c. 

.\ow'  that  the  governing  equation  system  for  the  wet  surface  unknowns  has  been  fully 
defined,  equations  (1)  are  rewritten  by  solving  the  first  for  x  and  substituting  into  the 
second.  In  combination  with  the  first  of  (1),  these  become 


(E//  -  E/,E,;E,/)p5  -  gy  -  E/,Ejg, 

-  E,/P5. 


(10) 


P5  is  found  from  the  first  of  (10)  while  x  is  then  obtained  from  the  second.  This  is  the 
solution  procedure  currently  implemented  in  the  SWEEPS  processor. 
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Section  3 

Governing  Equations  For  Far  Field  Pressure  Solution 


The  pressure  p  in  the  acoustic  field  is  given  by  the  Helmholtz  wave  equation 

V^p+k^p-^0,  (11) 

where  V“  is  the  Laplacian  operator.  A  convenient  boundary  integral  solution  to  (11) 
can  be  given  by  particularizing  KirchholT’s  retarded  potential  formulation  [12,13]  to 
the  steady  state  case  to  obtain  the  Helmholtz  integral  equation 

Psca  ^  ^2  ^  ikr){p,  -  ps)  e  (12) 

'iTT  J  [  r  on 

where  the  explicit  time  dependence  has  been  omitted,  and,  ps^ca  is  the  scattered  (or, 
radiated  pressure  if  pi  0)  at  any  point  P  outside  of  the  wet  surface  boundary  B,  r 
is  the  distance  from  P  to  a  point  Q  on  B,  n  is  the  outward  unit  vector  normal  to  B 
at  Q.  and  dr’dn  is  the  cosine  of  the  angle  between  r  and  n.  It  should  be  emphasized 
that  (12)  is  an  (  xact  result  and  that  the  only  approximation  involved  here  is  due  to  the 
fact  that  T  and  ps  are  obtained  on  the  wet  surface  by  the  DAA2  as  described  in  the 
proceeding  section.  A  geometric  approximation  is  also  invoked  when  (12)  is  discretized 
by  assuming  that  x.  P],  and  pg  are  constant  over  each  boundary  element  covering  the 
surface  B.  The  result  becomes 


.  A'  ,  r  1  . 

Psca  =  -  +  ^  (1  +  ikRi){pii  +  Psx)  e' 

'  1^1  '  L  * 


where  A'  is  the  total  number  of  boundary  elements  on  the  wet  surface.  A,  is  the  area  of 
the  element,  while  5,  and  'y,  have  already  been  defined  in  Section  2,  except  that 
here  they  pertain  to  the  point  at  which  the  scattered  pressure  is  to  be  calculated,  rather 
than  the  incident  wave  source.  (13)  then  forms  the  basis  of  the  T.\RGET  processor. 
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Section  4 

Submerged  Spherical  Shell  Solutions 

The  nondimensional  equations  of  motion  for  the  axisymmetric  modal  vibrations  of  a 
submerged,  linearly  elastic  spherical  shell  are  briefly  summarized  here  from  [3].  With 
the  addition  of  the  right  hand  side  forcing  terms  due  to  an  incident  plane  wave,  they 
may  be  expressed  in  matrix  form  as 


'An  -  n(n  r  l)tZ'^ 

Bn 

0  ■ 

'  Vn  ' 

0 

Bn 

Cn  - 

< 

Wn 

>  =  < 

-  Pn) 

0 

Rn{<^)- 

Pn  . 

in  which  v„  and  tu„  are  the  meridional  and  radial  components  of  modal  shell  displace¬ 
ment,  respectively,  p^,  p^  and  p^  are  the  scattered,  internally  forced  and  incident 
components  of  modal  surface  pressure,  respectively,  and  is  the  fluid  particle  velocity 
due  to  the  incident  wave,  p  —  pajp^h,  where  p^  is  the  density  of  the  shell  material, 
and  a  and  h  are  the  shell’s  radius  and  thickness,  respectively.  In  addition 

An  =  7^(1  +  «)”(”  +  ^)sn, 

+  tCn)n(n -L  1),  (15) 

Cn  =  'r^!2(l  +u)  +  n{n  -r 


Here,  q  -  c^/c.  c  =  h^/l2a^.  and  =  «(n  -h  1)  -  (1  -  u),  where  and  i/  are  the  plate 
velocity  and  Poisson's  ratio  for  the  shell  material,  respectively,  and  n  is  the  modal  index 
for  meridional  expansion  in  Legendre  polynomials  of  the  shell  displacement,  external 
pressure,  and  internal  pressure  fields. 

The  incident  pressure  and  particle  velocity  can  be  written  as  14 

=  Po(-l)"(2n  -  l)j„(i’). 

po  is  the  amplitude  of  the  incident  pressure  and  the  j„  are  the  spherical  Bessel  functions 
of  the  first  kind.  The  factor  (-1)"  occurs  because  the  incident  wave  is  a‘'''umed  to 
impinge  on  the  sphere  ai  0  =  n. 

For  the  D.\A2  solutions  required  here,  the  polynomials  (?„(!)  and  taiv  the 

following  forms  :3] 

DAA^  :  Qn(^)  =  *  fif(”  ^  !)■  +  7^'(n  -  1)  •  g(ti  ■  I)'. 

(17) 

DAA2  :  Qn(w)  =  iu)  n,  "  iSj{n  ^  1)  n(n  •  1). 
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Finally,  it  should  be  noted  that  the  nondimensionalization  of  (14)  is  based  on  the 
relations 


v  =  V/a,  w  —  Wja,  u  =  U/c,  p  =  P  f  pc^ ,  Cj  —  ujajc,  (18) 

where  V ,  IF.  U,  P  and  .t’  are  the  appropriate  dimensional  variables. 

The  wet  surface  solution  to  the  3x3  matrix  system  in  (14)  is  easily  obtained  numerically 
as  a  function  of  the  mode  number  n.  The  pressure  in  the  fluid  field  can  then  be  found 
by  using  this  solution  in  the  infinite  series  expansion 


p.ca{r,0)  -  i'Sj‘^  Pn{cos0)\{pi,  +  P^)i'(il')  -  i^u’nin(^')'(>n(i^r)  iyn(iZ'r)),  (19) 


wliere  r  and  0  are  the  radial  and  meridional  polar  spherical  coordinates,  respectively, 
and  the  nondimensional  radius  r  is  given  in  terms  of  the  dimensional  radius  P  as 
r  R  a.  Pn{cos6)  are  the  axisymmetric  Legendre  polynomials  and  y„  is  the  spherical 
Bessel  function  of  the  second  kind.  The  minus  sign  multiplying  yn  differs  from  jl4] 
because  the  time  dependence  used  here  is  e^^'^  rather  than  that  is  used  there. 

(19)  follows  directly  from  (4.14)  and  (8.10)  of  :14  ,  although  it  is  not  explicitly  derived 
there,  since  the  exact  solution  for  the  sphere  can  also  be  expressed  more  simply  in 
terms  of  the  modal  surface  displacements  or  pressures  alone.  This  is  in  contrast  to  that 
for  an  app)roxiniate  wet  surface  solution  like  the  DAA2  . 
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Section  5 

Code  Verification  Using  The  Spherical  Shell  Problem 

Two  problems  involving  analytical  solutions  for  the  subnu-rged  spherical  shell  are  de¬ 
scribed  in  this  section  for  the  purpose  of  providing  a  basis  of  verification  of  the  SWKKPS 
and  TARCKT  proc  essors.  The  first  involves  comparisons  of  the  internally  forced  modal 
results  of  3  with  those  produced  by  SWEKPS;  while  the  second  is  a  test  of  both 
S\\  KKPS  and  T.ARCET  for  an  incident  plane  wave  excitation  that  was  motivated  by 
ir>i. 

The  computational  rnodc'l  chosen  for  the  comparisons  to  be  made  here  is  a  cpiarter  of  a 
sphere  \shose  surface'  mesh  consists  of  cpiadrilatc'rals  only,  all  of  which  are  very  nearly 
scpiarc'.  It  cemtains  OG  c'lc'iiients  and  1  l.'t  node  points  and  is  sliown  in  Figure  1  in  an 
exf>lodc'd  view.  It  was  construe  tc'cl  using  a  pre-procc'ssor  dc'cc'loped  especially  for  this 
purpose  16  and  forms  the  basis  of  the  structural  clc  inc'iit  grid  for  ST.\GS,  and  also 
for  the  fluid  element  grid  for  the  I  SA  procc'ssor  FIT  M.\S. 

Tiiferiinlly  Forced  Problem 

In  this  case  the  analytical  DAAf  calculations  of  3  are  used  to  provide  the  basis  for 
comparison  with  SWEEPS  predictions  of  surface  radial  \eloc  ities  and  surface  pressures. 
The  modal  \ibration  results  presented  in  .3  include  plots  of  these  quantities  as  func¬ 
tions  of  frequency  and  clearly  show  the  resonance  and  anti-resonance  zones  that  are 
of  primary  interest  in  this  probloiri.  The  particular  computations  reported  there  are 
carried  out  for  the  modal  indices  n  ■  0, 1,2  and  3  and  for  two  non-zero  values  of  struc¬ 
tural  daiiiping.  In  contrast,  SWEh^PS  does  not  currently  include  structural  damping 
and.  of  course,  is  constructed  around  a  discrete-element  model  that  can  include  many 
mocles.  'I'o  make  the  comparisons  then,  the  governing  equation  syste-n  presented  in 
3  and  stirnmari/ed  in  the  preceeding  section  is  solved  for  a  selected  number  of  fre- 
(|uencies  without  structural  damping.  Also,  the  surface  distribution  of  the  SWEEPS 
internal  forcing  function  is  specified  to  be  that  of  the  appropriate  axisymmetric  Legen¬ 
dre  polynomial  mode.  The  results  presented  in  Figures  2  through  9  correspond  to  the 
parameter  values  used  in  3 

a  h  100.  p,.., V  7.67,  u  0.3.  r.,  c  3.53. 

and.  as  can  be  seen,  the  comparisons  are  excellent  except  for  one  point  at  u;  =  3  for 
n  0.  .No  immediate  explanation  for  this  minor  discrepancy  can  be  offered  at  this 
t  irne. 

luciclcmt  Plane  Wave  Excitation 

The  work  reported  upon  in  15  is  a  feasibility  study  that  examines  the  accuracy  of 
several  surface  interaction  approximations  as  applied  to  the  underwater  acoustic  echo 
signal  problem.  7'hat  work  makes  use  of  the  known  exact  and  approximate  modal 
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Figure  10  Backscattered  Pressure  At  r  =  20  For  Submerged  Spherical  Shell 
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Section  6 

Discussion  and  Conclusions 

Based  upon  the  excellent  agreement  between  the  discrete-element  and  closed  form 
DAA2  solutions  obtained  here,  it  may  be  concluded  that  the  SWEEPS  and  TARGET 
processors  are  capable  of  reproducing  the  variety  of  phenomena  that  are  inherent  in 
underwater  acoustics  problems.  In  addition,  since  it  has  already  been  demonstrated 
elswhere  that  the  DAA2  approach  provides  very  good  approximations  to  exact  solu¬ 
tions,  it  is  apparent  that  such  discrete-element  formulations  will  be  valuable  assets  in 
the  acoustic  studies  of  submerged  structures. 

The  results  of  this  paper  do  raise  an  interesting  question  in  that  the  discrete-element 
computations  seem  to  retain  accuracy  to  higher  frequencies  than  would  be  expected. 
The  reason  for  this  serendipitous  behavior  is  not  yet  understood. 

With  regard  to  software  details,  SWEEPS  and  TARGET  currently  have  some  minor 
restrictions  on  their  usage.  Common  to  both  is  the  requirement  of  an  infinite  fluid, 
ie.,  the  presence  of  a  free  surface  is  not  treated.  In  addition,  structural  damping  has 
not  yet  been  implemented  in  SWEEPS.  The  inclusion  of  these  two  capabilities  requires 
some  additional  but  straightforward  code  enhancement. 

.At  this  time,  the  computational  algorithm  for  SWEEPS  is  a  direct  elimination  solution 
of  the  structural  and  iluid  equations,  hence  the  CPU  time  is  a  linear  function  of  the 
number  of  frequencies  desired.  It  is  also  roughly  proportional  to  N^,  where  N  is  the 
number  of  fluid  degrees  of  freedom.  Sample  execution  times  for  a  single  frequency  for 
57,  113,  and  270  fluid  DOF  problems  are  1,  5,  and  60  CPU  minutes  respectively  on 
the  V.AX  11/780.  In  order  to  reduce  such  expense  for  multi-frequency  calculations, 
an  incremental  iterative  scheme  is  planned  that  can  eliminate  costly  complex  matrix 
factorizations  at  every  frequency.  It  should  also  be  emphasized  that  an  increase  in 
efficiency  in  the  SWEEPS  processor,  through  a  change  in  the  algorithm,  is  possible  only 
because  of  the  frequency  independent  matrices  that  are  inherent  in  the  DA  A  approach. 
Such  an  algorithm  would  have  no  effect  upon  current  exact  solution  methods. 
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solutions  for  the  spherical  shell  irradiated  with  a  plane  incident  wave  similar  to  those 
given  in  the  preceeding  section  (see  [3]  for  the  exact  polynomials  and  P))- 

In  particular,  the  backscattered  pressure  is  computed  at  r  =  20,  6  -  n  and  is  plotted 
as  a  function  of  frequency  in  terms  of  the  echo  function  defined  as 

(20) 

Po 

It  concludes  that  the  DAA^  with  g  -  1  is  a  very  good  approximation  to  the  exact 
solution  over  the  entire  frequency  range  studied  there:  0  <  <  16. 

For  the  purposes  of  the  current  work,  comparisons  will  be  drawn  between  discrete- 
element  and  closed  form  DAA™  solutions  using  (19)  of  the  preceeding  section  rather 
than  make  use  of  any  results  already  presented  in  15'.  The  reason  for  this  is  that  the 
scattered  pressure  in  the  fluid  was  obtained  only  aj)  roximately  there,  using  numerical 
integration  of  the  Helmholtz  integral,  and  some  differences  were  clearly  evident  between 
those  and  the  current  discrete-element  results,  particularly  in  the  range  2  <  u.’  <  6. 
Using  (19)  then,  the  P„(coii6)  reduce  to  (  1)"  for  the  case  0  -  tt  \n  order  to  compute 
the  backscattered  pressure.  In  addition,  the  convergence  criterion  is  chosen  so  that  the 
last  modal  increment  to  ;/oc!  is  less  than  10  “  times  the  current  value  of  j/ool-  The 
results  presented  in  Figure  10  correspond  to  the  parameter  values 

a,'h  -  39.5,  pg/p  —  2.7,  t/  0.355,  c^/c  =  3.7928, 

that  were  used  in  |15;.  As  can  be  seen,  the  agreement  is  excellent  over  most  of  the 
frequency  range  and  it  is  only  when  the  surface  discretization  is  very  coarse  in  relation 
to  the  wavelength  of  the  incident  pressure  that  the  SWEEPS  calculations  begin  to 
diverge  wildly  for  >  14.  Indeed,  it  is  surprising  the  results  are  so  good  even  out  to 
w  =  12,  since  at  this  frequency  there  are  only  3  mesh  points  per  wavelength. 

It  is  easily  demonstrated  that  these  results  are  not  particularly  dependent  upon  well 
defined  low  order  modes  that  have  a  sufficiently  large  number  of  mesh  points  on  a 
modal  wavelength.  The  convergence  of  !/oc!  is  shown  in  Figure  1 1  as  a  function  of 
mode  number,  and  it  is  seen  that  n  —  12  is  still  a  significant  contributor.  For  the 
discrete-element  solution,  this  mode  has  slightly  less  than  3  points  per  wavelength 
around  the  sphere. 

Incidently,  since  the  results  of  Figure  10  agree  so  well  for  two  so  completely  different 
numerical  methods,  it  would  appear  that  there  is  a  minor  numerical  flaw  in  some  of 
the  results  of  jlSj. 
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